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Abstract 

We study the emissivity properties of a geometrically thin, optically thick, steady accretion disc 
about a static boson star. Starting from a numerical computation of the metric potentials and the 
rotational velocities of the particles in the vicinity of the compact object, we obtain the power per unit 
area, the temperature of the disc, and the spectrum of the emitted radiation. In order to see if different 
central objects could be actually distinguished, all these results are compared with the case of a central 
Schwarzschild black hole of equal mass. We considered different situations both for the boson star, 
assumed with and without self-interactions, and the disc, whose internal commencement can be closer 
to the center than in the black hole case. We finally make some considerations about the Eddington 
luminosity, which becomes radially dependent for a transparent object. We found that, particularly at 
high energies, differences in the emitted spectrum are notorious. Reasons for that are discussed. 
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P-il Introduction 

Oh: 

^ The possible existence of very massive non-baryonic objects in the center of some galaxies is being studied 
since a long time. As far as we are now aware, they were first hypothesized by Tkachev Q, who also 
studied the emissivity properties arising from particle anti-particle annihilation processes. More recently, 
^ detailed studies of the properties of neutrino ball scenarios have also been carried out by Viollier and his 

^ collaborators |||]. In Ref. in addition, we have also explored whether supermassive non-baryonic boson 
■ " " stars might be the central object of some galaxies. To fix the situation to a particular case, we have paid 
special attention to the Milky Way. This study had a twofold aim. On one hand, it focused on what current 
dynamical observational data have established regarding the properties of the galactic center. We have 
concluded in this sense that scalar stars fitted very well into these dynamical constraints. On the other 
hand, we have also discussed what kind of observations could actually distinguish between a supermassive 
black hole and a boson star of equal mass, pointing out several possible tests. 



In the case of our own Galaxy, recent observations probe the gravitational potential at a radius larger 
than ~ 10^ Schwarzschild radii [Q]. The black hole scenario is the current paradigm, but suggestions 
that the dark central objects are indeed black holes are based only on indirect astrophysical arguments, 
basically dynamical in nature, that could be sustained by any small relativistic object other than a black 
hole, if it exists ||5|. It is then advisable to explore possible alternative scenarios. The final aim should be 
to devise definitive tests that could observationally solve the issue. One is to look for the event horizon 
projected onto the sky plane. Although this is the key concept of both, Constellation-X |^] and Maxim |^ 
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satellites, the needed spatial resolution is still a dream for the future. It might be more feasible to look for 
the realization of the recently introduced concept of the shadow of a black hole § . Although any highly 
relativistic object would also produce a shadow, the observed features might be enough to decide whether 
an event horizon is present or not. 



Prom the point of view of boson stars physics (see for a recent review), we also need to know whether 
fundamental scalars capable to form the stars do exist; observational tests in this sense are highly desirable 
too. Only after the discovery of the boson mass spectrum we shall be in position to determine which 
galaxies, if any at all, could be modeled with such a center. In recent works, Schunck and Liddle ||lO|| , 
Schunck and Torres |11], and Capozziello et al. |12|, among others, analyzed different observational effects 
that boson stars would produce. Also, boson stars were proposed as sources for some of the gamma ray 
bursts 1 13 1, and as a possible lens in a gravitational lensing configuration |14], following recent interest in 
analyzing the gravitational lensing phenomenon in strong field regimes [15]. 



However, as far as we know, literature lacks the study of an accretion process, onto a boson star, either 
massive or supermassive. Do the emissivity properties of the disc differ when the central object, instead 
of being a black hole, is assumed to be a boson star? Can we detect these differences? From where these 
deviations, if any, come from? Can accretion onto boson stars help model galactic centers? Of what kind? 
To completely answer these questions would require the analysis of different models of accretion discs, with 
various degrees of complexities. In this first approach, we shall analyze the properties of the simplest, 
steady, geometrically thin, and optically thick accretion disc model, rotating onto a static boson star. We 
shall compare all our results with those obtained using a black hole of the same mass. 

The fact that all circular orbits are stable for static stars (as we shall show), can pose a problem to 
accretion scenarios upon non-baryonic (with no surface) static objects. Accretion would follow a series of 
stable circular orbits, loosing angular momentum and radiating part of the generated heat. If particles 
can always found a stable orbit, provided an enough amount of time, they would all end up in the center, 
and should a way of diverting them from there not exist, we would confront the formation of a baryonic 
black hole in the center of every non-baryonic star subject to overdense environments. This problem have 
apparently (as far as we are aware [|) been not clearly mentioned ever before in the literature of boson 
star solutions (see [^), although we have no other option than confront it if we are to talk about any 
astrophysically relevant use of these scalar star models. 

This problem can, however, be alleviated in more general situations. In the rotating case, particles 
orbits were analyzed by Ryan |jl^ (see his section IV). He has shown that circular geodesic orbits are not 
stable beyond a given point, located at about 5/3 times the radius of the doughnut hole which appears 
in rotating boson star solutions. He has considered the swirling of a stellar size object (a black hole, or 
neutron star) within a supermassive rotating boson star and studied the gravitational wave emission as 
a mechanism for detection. Accretion would not continue in circular orbits since there is none after that 
point. Also, we are just analyzing the case in which a particle continues to travel in geodesies within the 
star interior, disregarding any possible influence of the boson star matter. And there is also the fact that 
two body encounters will be unavoidable in the innermost regions of the boson star, due to the increased 
matter density. These two body encounters will shift the particles to superior orbits where they'll find 
turning points, and bounce. 

But even if an inner black hole forms, the influence upon the accreting matter that it would exert would 
^We acknowledge e-mail discussions with Dr. Tkachev in early 2000 on issues related to this point. 
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be limited by the amount of mass it has compared with the mass of the non-baryonic object (also note that 
the boson star radius and the shells where most of the non-baryonic matter is located, are farther away 
than the Schwarzschild radius of the presumed black hole by a factor of at least 100). In general, there 
will be situations in which the accretion rate will be so low, that even if a black hole is formed with all the 
mass accreted during the lifetime of the universe, it will still have several orders of magnitude less than 
the boson star mass. Consider the center of the galaxy, its mass is believed to be above 2 10^ Mq, while 
the accretion rate is ~ 1O~^M0 yr~^. Then, in the absolutely worse case, if a black hole is formed with the 
accreted mass in a period of lO^'^ yr, it will have 10^ Mq, and its gravitational influence will be defied by 
the non-baryonic object. Boson stars containing fermion objects within have been considered in the past 
|17|, and this appear to be an extreme case where the fermion (neutrons) component have collapsed to a 



black hole. 

Finally, a detailed analysis of the evolution of boson stars subject to continuous inflow of non-baryonic 



particles was carried out in Ref. |18|. Their results showed that under finite perturbations, the stars on 
the stable branch will settle down into a new configuration with less mass and a larger radius. Then, 
the accretion of non-baryonic matter possibly entering into the condensate would not pose a problem to 
boson star stability, nor generate a collapse. We recall that the instability of a boson star with respect 
to gravitational collapse has been studied by Kusmartsev et al. (see Ref. see also Ref. [|^]), using 



catastrophe theory. Rotating boson stars were studied (among others) by Schunck et al. |21 



In summary, even when the physical model used here could be regarded as too simplified and not 
complete (the star is not rotating, and a mechanism by which diverting matter from the center is not 
specifically given) we believe it still is important to begin to address the issue of real astrophysical scenarios, 
as accretion, upon theoretically foreseen non-baryonic objects. One of the first steps in this direction is 
presented in this work. 



2 Steady and geometrically thin accretion disc 

A test particle rotating around an spherically symmetric object, which generates a space-time described 
by the metric 

ds^ = e'^^^'Ut^ - e^('^)dr2 - r'^{dd'^ + sin^ d d<f^) , (1) 

would do so at a velocity v'^ = rv' jl [^.0 A prime stands for the derivative with respect to the radial 
coordinate. Let us call, for the ease of notation, B = e'^ , then B' / B = v' . In the case of a black hole, B is 
given by the Schwarzschild metric B{r) = 1 — 2M/r, what finally yields the Keplerian velocity = M/r. 
The angular velocity, defined as usual as Vl = v^/r, is then proportional to r~^/^ and it decreases rapidly 
far from the center. For a more generic case, the rotational velocity and its derivative will be given by 

(2) 



17 = W—, Vt' 



We now turn to dimensionless quantities. The radial coordinate being x = mr, where m will be 
associated with the mass of the boson. Then, v"^ = dB/dx x/2. Recall that c is being implicitly taken 
equal to 1, the dimensionfull velocity will just be w,^ x c. Equivalently, the dimensionfull rotational velocity 

^This form of the velocity (and the metric) may not be generally valid in the region dubbed as "dark matter zone" in the 
paper by Nucamendi et al. fed], i.e. where the flat rotation curves of galaxies are found. 






-1/2 


B" 


B'- 






.27 


~ 2^. 



3 



will he Q = m c y^[{dB/dx)x/2], and 



Ub(x) X ,,a.5228 X 1024 
^2"^[^^^] ^ ■ 



n' 



dn 2 7.72973 x lO^^ 



dx 



m^[GeV]- 



km s 



(3) 



The power per unit area in the black body disc we are studying is given by [34| 
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where Ri is the interior limit of the disc and fij is the velocity there. M is the accretion rate, it is an 
external parameter that should be observationally obtained, Ri should also be either observationally derived 
or theoretically estimated, by taking into account the nature of the central object, all other quantities are 
known if the explicit form for the metric coefficients is. As we see from Eq. (^), the standard model 
for a steady, geometrically thin and optically thick disc predicts a surface emissivity that is independent 
of viscosity, given all other disc parameters. In fact, Eq. (^) is strictly true only when the gravitational 
potential is ~ GM/r, otherwise relativistic corrections would appear. The benchmark case for disc accretion 
is that of a Schwarzschild black hole. For it, M is independent of r, it is just a constant equal to the black 
hole mass, Mbh, and it is easy to analytically derive the form for the rotational velocity and its derivative 
by using the expressions given above. With these functions, the power per unit area for a central black 
hole of mass Mbh is 



ovr 



1- ( ^ 
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x-^M[MQ/yr] m'^[GeY] 1.47065 x 10 



74 eJ^g 
cm^s 



(5) 



The appearance of m in the latter formula is justified because of the use of the dimensionless radial 
coordinate and the definition of MBH,a by means of the equality Mbh = MBH,a Mpjm, where Mpi is 
the Planck mass. This particular form of writing the hole mass is useful because it will allow a direct 
comparison with the case of a boson star. In a general situation (when not necessarily a black hole is the 
central object), the power per unit area will be given by 



D{r) 



BS 



1 



1 



xM[MQ/yr] m?[Ge\] 1.47065 x 10 



74 erg 



cm^s 



(6) 



Here, [il] and [0'] stand for the dimensionless form of these quantities. The superscript BS refers to our 
later use of this formula for the case of a boson star. 



If the disc is optically thick, as we are assuming, the local effective temperature must be sufficient to 
radiate away the local energy production. The temperature of the disc is then related with the latter result 
by D{r) = crT^, as it is appropriate for a black body system; a being the Stefan-Boltzmann constant (equal 
to 5.67 x 10~^ erg s~^ cm"^ K"^). For the luminosity, L{v), and flux, F{v), we shall use the relationship 

L{u)=And'Fiu) = J^^ ^../.T _ 1 > (7) 

where d is the distance, Rout is the outer border of the disc, i the disc inclination, and h{= 6.6256 x 10"^"^ 
erg s) and k{= 1.3805 x 10~^^ erg K~^) are the Planck and Boltzmann constants, respectively. In order 
to make the first comparison we shall take Ri = 3Rg = 3 x 2GM = QM~i^ M^i^ = 6Ma, and 

R^^ = 30 - 50Ep = 30 - 50 X 2GM = 60 - lOOAf-^M^, ^ ^ Xout = 60 - 100M„, with M„ being the 
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dimensionless mass. This will let us know whether, even in the non-relativistic regime and for the same 
initial commencement of the disc, we are able to see any difference in the emitted spectrum. The relativistic 
analysis will follow in the next section. 



In dimensionfull form, for an inclination of 60 degrees, changes of variables yield a luminosity equal to, 
= 2.2559 X 10-^3,rg rn-^[GeV] .3[,-3] ^^^^V^tR^ ■ (») 

^hu/kT{x/m) ^y^^^ obtained for each particular central object, through the solution for the metric coeffi- 
cients and the computation of the rotational velocities and temperatures (for the latter we use D{r) = crT^, 
with the corresponding power) We shall now apply the presented formulae to the case of a boson star. 
For readers not familiar with the way in which boson star configurations are obtained, we suggest Refs. 



24, H, 27, 19, 25, 29, BG] 



3 Numerical results for a non-relativistic disc 

We can now note the following basic steps of the computation of the emission properties of the disc: 1) 
Define the boson mass m, the inner and outer borders of the disc and the accretion rate M; 2) Compute 
^hu/kT(x/m) £qj, (jjffgj-gn^ frequencies; 3) Integrate /^f™' ^hiy/kT{x/m)_i obtain the luminosity for a given 
frequency; 4) Repeat the previous steps for several frequencies and build up the spectrum. We shall take 
a model with central mass equal to 2.5 x W^AIq. ^ We shall take M to be equal to a few x lO~^M0yr~-^; 
in the computation that follows, if not otherwise stated, the pre-factor will be taken as 2. In this situation, 
Xi will be equal in this case to 3.798 and Xout = (30 — 50) Xg, with Xg = 1.266. 

Remark on scaling: Note that if we would be interested in other different masses for the compact object, 
we would just need to change the boson mass. All quantities are already scaled, so that a change in m, 
and eventually in M, will readily give the values of all other physical results. 

In Figure |l| we show the dimensionless rotational velocities, both and fi, for the previously quoted 
model. Through the comparison of the rotational velocities that a black hole of equal mass would produce, 
it is easy to see that a boson star is a highly relativistic object. The velocities for a black hole center are 
also shown in the Figure. A comparison of the power per unit area is shown in Figure |2[ It is also shown 
there the difference between the temperatures of the disc. As one would expect from the behavior of the 
metric coefficients, the boson star curves, especially farther from the central object, tend to mimic those 
of the black hole case. However, in the inner parts of the disc, it is apparent that a slight deviation can be 
noticed. The boson star produce more power per unit area, and a hotter disc, than a black hole of equal 
mass. In order to see whether the dependence on the self-interaction parameter is strong, we have computed 
these same quantities for the case of A = 100 and the same total mass and accretion rate. We have used a 
central density equal to cr(0) = 0.093, what yields a boson star mass -in dimensionless units- equal to 2.25. 

■^Since we have already assumed that A = 0, for a maximal boson star with central density (t(0) = 0.271, we shall need the 
boson mass to be equal to ~ 10"^'' GeV. A discussion of this value of m on the light of particle physics was given in Ref. 
see also Ref. In any case, note that because of the scaling properties, the general trend of the results presented here is 
valid for every value of m, even for those leading to stellar sized boson stars. This ~ 10~^^ GeV value is obtained equating 
the boson star mass, Mbs = 0.633 M^Jm to the mass of the central object, where, in the previous equation, 0.633 = MBS,a 
is the value of the boson star mass in dimensionless units; this value is numerically obtained solving the star equations of 
structure. In Ref. [^, we presented a relation between the boson mass m expressed in GeV and the total mass of the star 
expressed in millions of solar masses. This relationship is m [GeV] ~ 1.33 X 10-'^^MBS.a/M[10'^MQ], and was used to obtain 
the value of m quoted above. 
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Figure 1: Left: Dimensionless velocities. The upper set of curves are v^f,, the lower one are the rotational 
velocity 0, both for a boson star (solid) and a black hole (dash). The boson star has a dimensionless 
radius equal to xbs = 4.3, and a dimensionless mass equal to 0.633. The black hole has equal mass. Right: 
Dimensionless derivative of the rotational velocity, $7'. 

In order for this dimensionless mass to rightfully represents the size of the central object, we need a boson 
mass equal to 1.20 x 10~^^GeV. The final spectrum of the radiation emitted by the disc was computed 
and it is shown in Figure |3|. It is clearly seen that throughout most parts of the electromagnetic spectrum, 
the obtained differences in the disc properties, such as the temperature or the power per unit area, will 
not be observationally noticed. Differences between the boson star cases with and without self-interaction 
are negligible. However, it appears that differences between black hole and boson stars centers turn out 
to be important at the most energetic part of the electromagnetic range, beginning in the far ultraviolet. 
The disc model we are considering does not produce energy in the X-ray and gamma-ray bands, but the 
possibility is open that even in the simplest non-relativistic case we have considered, a boson star and a 
black hole case could in principle be distinguished by the radiation properties at high energies. We can 
expect that better and more realistic models of accretion discs can show stronger differences. 

The metric potentials of a boson star make the rotational velocities and their derivatives, except for 
the innermost parts of the disc, always greater than their counterparts produced by a central black hole. 
Looking at the left panel of Figure ||, it can be seen that only for the more central coordinates, the black 
hole rotational velocity is greater than that produced by a boson star. From the right panel we see, in 
addition, that the derivative of O is also always greater for a boson star than for a black hole. Recalling 
that D{r) cx ilil', a greater rotational velocity and rotational velocity gradient implies a greater power. 
The same happens for the temperature of the disc. The innermost part of Figure |l], where the rotational 
velocity of the particles rotating around a black hole diverges, can not be taken into account. These par- 
ticles are all following unstable orbits and are not part of the disc, which begins at Xi = 'ixg, at a value of 
about 4 in the x-axis of Figure |l|. The effect of a larger temperature and power per unit area is then the 
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Figure 2: Left: Dimensionfull power per unit area for a black hole, and boson stars with self-interaction 
equal to and to 100, all of the same mass. The position where the maximum of the curves are is zoomed 
out in the figure. Right: Dimensionfull temperatures of the disc for a black hole, and boson stars with 
self-interaction equal to and to 100, all of the same mass. The position where the maximum of the curves 
are is zoomed out in the figure. 

integrated output of a disc rotating faster when a boson star is the central object. 

Before ending this section we shall provide a brief discussion on the boson mass chosen. This discussion 
is based in our previous paper |^], as well as in 1 33]. As discussed in Q, based on the constraints imposed 
by the mass-radius relationship valid for the scalar stars analyzed, we may conclude that 1. if the boson 
mass is comparable to the expected Higgs mass (hundreds of GeV), then the center of the galaxy could be 
a non-top ological soliton star p6| ], 2. an intermediate mass boson could produce a super-heavy object in 
the form of a boson star, and 3. for mini-boson stars to be used as central objects for galaxies the existence 

^We have been using central object and accretion disc parameters, like mass and accretion rate, consistent with those 
inferred for our own Galaxy We can then ask if the replacement of the presumed central black hole for a boson star can 
do any better in fitting the observationally obtained spectrum. The answer to this question is no. However, this is not (at 
this time) because of the fact of the different central object, but because of the properties of the accretion disc itself. The disc 
model here considered simply does not work for Sgr A*, neither with a black hole, nor with a boson star. The measured M 
and M for Sgr A* would yield -within the standard theory of a steady, optically thick, geometrically thin accretion disc- an 
accretion luminosity equal to O.lMc^ > 10*" erg s~^, assuming a nominal efficiency of 10%. However the total luminosity of 
Sgr A* is less than 10"^^ erg s~^. In addition, the value of M would make the standard accretion disc broad band spectrum to 
have its peak in the infrared region, what is opposite to observation. The spectrum is -with the exception of a few bumps- 
essentially flat, with detected flux even in x and 7-rays. A recent compilation of luminosity measurements for Sgr A* was given 
by Narayan et al. [^. Observations suggest, then, that Sgr A* is not behaving like a blackbody, disregarding the nature of 
the central object. An alternative solution for the blackness problem of Sgr A* is the so-called advection dominated accretion 
discs, or ADAF | ^ , ^ . The key for ADAF models is that most of the power generated by disc viscosity is advected into the 
hole, while only a small fraction of it is radiated away. An essential ingredient for ADAF models is then the existence of an 
event horizon. We remark that for the less active galaxies, like ours, a boson star seems to be a more problematic center than 
a black hole: the blackness problem could even be more severe. 
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Figure 3: The final emission spectrum. Differences between tiie boson star cases witfi and witfiout self- 
interaction are negligible, they are not noticed in the plot. Noticeable deviations between a supermassive 
black hole and a boson star begin at the far ultraviolet. 

of an ultralight boson is needed. These conclusions are to considered as order of magnitude estimates. 
If boson stars really exist, they could be the remnants of first-order gravitational phase transitions and 
their mass should be ruled by the epoch when bosons decoupled from the cosmological background. The 
Higgs particle could be a natural candidate as constituent of a boson condensation if the phase transition 
occurred in early epochs. A boson condensation should be considered as a sort of topological defect relic. 
If soft phase-transitions took place during cosmological evolution e.g., soft inflationary events, the leading 
particles could have been intermediate mass bosons and so our super-massive objects should be genuine 
boson stars. If the phase transitions are very recent, the ultralight bosons could belong to the Goldstone 
sector giving rise to miniboson stars. We should also mention the possible dilatons appearing in low-energy 
unified theories, where the tensor field of gravity is accompanied by one or several scalar fields. In string 
effective supergravity, the mass of the dilaton can be related to the supersymmetry breaking scale msusY 
by ~ lQ~^{msusY /TeV)"^ eV. Finally, a scalar with a long history as a dark matter candidate is the 
axion, which has an expected light mass m„ ~ 7.4(10'^ GeV//^)eV ~ 10"" eV with decay constant f„ 
close to the Planck mass. 

4 Coordinate dependent Eddington luminosity 

It is interesting to make now a comment on the Eddington luminosity {Le) when a boson star is the 
central object. Consider first the basic concept. If the luminosity produced by the accreted material is 
too great, then the radiation pressure it would produce would blow up the infalling matter. The limiting 
luminosity, known as Le, is found by balancing the inward force of gravity with the outward pressure of 
radiation [^]. This limit is then found by assuming that the infalling matter is fully ionized, and that 
the radiation pressure is provided by Thompson scattering of photons against electrons -which in turn are 



8 



E 1E42 




■ Eddington luminosity 
Disc luminosity HAC 
Disc luminosity LAC 



5 10 15 

Gravitational radii (2GM) 



20 



Figure 4: Comparison between the dimensionfull Eddington luminosity, now dependent on the radial 
coordinate, and the power radiated by the disc at different radius. HAC stands for a high accretion rate, 
here taken to be equal to 2 x lO^^MQyr^"*^. LAC stands for a low accretion rate, equal to 2 x lO^^MQyr""*^. 
The coordinate dependent Eddington luminosity is everywhere greater than the power generated by a thin 
disc about a boson star. Note that the form of the curve for the Eddington luminosity is provided by that 
of the mass distribution within the boson star. 



strongly coupled with the protons resident in the plasma by the electromagnetic interaction. Equating 
the involved forces results, in a first approximation, aTL/Anr^c = GMrrip/r'^ , where ax is the Thompson 
scattering cross section and rrip is the proton mass. When this relationship is fulfilled, L = Le- Le is the 
maximum luminosity which a spherically symmetric source of mass M can emit in a steady state. Longair 
|l35|| provides the value of the Eddington luminosity in watts by introducing the gravitational radius Vg, 
Le = 2-Kr gVapC^ / ut = 1-3 x IO^^M/MqW. None of the active galactic nuclei were found to exceed this limit 



|35|. It is interesting to note that if M is constant, as in the case of a black hole, this limit is independent 
of the radius, that is, if L < L^; at any given radius, the inequality will be sustained everywhere. A boson 
star -as well as any other transparent object- has a non-constant mass distribution and the Eddington 
luminosity become a coordinate dependent magnitude.^ In this situation, one may ask if there is any point 
within the stellar structure such that, even if L < Le outside the star, the opposite inequality is valid at 
this point. If this is so, and L > Le^t = vq), steady accretion would not proceed. We would so be able 
to define the internal border of the disc. LE{r) is given by the distribution of mass, since LE{r) oc M{r). 
For the luminosity produced by the disc we have to integrate the power produced by a disc slab of size dr 
from Tj to r. This is what should be compared with LE^r). Then, we compute the total power radiated 



between r,- and r as 



/r Jut py„ rr 1 

D{r)2TTrdr = -5.705 x 10^^ M[— ^] ^ J 2 ^^'^ 



boson star could well turn into a non-transparent object if we are to admit the possibility of scalar electrodynamics 
effects, where considerable cross sections with photons may appear. This is an interesting problem, yet to be attacked, of 
boson star physics. 
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In Figure ^, we show the results for the Eddington luminosity for a non self-interacting boson star, compared 
with the power radiated by the disc for two different accretion rates. In none of these cases, the power 
generated in circles whose radii are smaller than the star size overcomes L^- This also supports the idea 
that accretion can continue inwards, within the boson star structure. 



5 A relativistic treatment for the accretion 



Even when useful as a first approach, the presented treatment fail when considering the innermost part 
of the disc, especially when the disc orbiting a boson star has a more internal commencement than that 
orbiting around a black hole. In order to get a complete picture, we have to treat the problem in a 
relativistic way, this is objective of the rest of this paper. 



5.1 Particle orbits 



Consider again the stationary, spherically symmetric, time independent, line element, written as 



-B{r)dt' + - 



1 



2M(r)/r 



dr^ + r^(d6'^ + sin^ 



(10) 



Again, the Schwarzschild solution is a special case of Eq. ([lO|) , occurring when B{r) = 1 — 2M{r)/r and 
M(r) is a constant equal to the mass of the black hole. In our present situation, B{r) must be consistent 
both with an asymptotically flat space-time and with the absence of a singularity. M(r), in turn, decreases 
with decreasing radius, being equal to at r = 0. The particle orbits in that general metric will be deter- 
mined by the conserved quantities, E = —pt : total energy; L = pg : component of the angular momentum. 
An additional constant of motion will be the mass of the test particle, let us call it m, which can be ab- 
sorbed by redefining quantities in a per-unit-mass basis. The general (non-trivial) equations for the orbital 
trajectory in a space-time described by Eq. (|^) are: 




nl/2 



(11) 



dX 



L 



sm 



(12) 



For a general derivation see the book by Weinberg [^] (take caution with the different definition of the 
symbols). Here, A is an affine parameter related with the proper time r, by A = r/m. We shall take m = 1 
for massive particles. A prime will denote derivation with respect to r. 



We are interested in the circular orbits. For such class, dr/dX must vanish instantaneously and at all 
subsequent moments, what imposes the conditions, V{r) = and V{ry = 0. These equations determine 
E and L, which per unit mass, are given by 

rh~ ~\2B-rB') ' m ~ ~ \2B - rB' ) ' ^ ^ 

These quantities are shown -both for the black hole and the boson star case- in the left panel of Figure 5. 
One can immediately show that both E^ and reduce themselves to the Schwarzschild expressions when 
B(r) = 1 — 2M/r, and that they tend to them asymptotically. Circular orbits might not exist for all values 
of r. It is needed that the denominator of E^ and L^, be well-defined, i.e. 2B — rB' > 0. But in the case of 
a boson star, this happens for all values of the radial coordinate, see the right panel of Figure 5. Then, in a 
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5 10 15 5 1015 20 25 

Dimensionless coordinate Dimensionless coordinate 

Figure 5: Left: Radial dependence of the energy, angular momentum, and the difference between them 
-with the momentum weighted by the rotational velocity fi- for circular geodesies, all per unit mass. The 
divergent curves are the corresponding magnitudes for the Schwarzschild black hole case. Right: Radial 
dependence of the denominator of the expressions for and L"^. Being the values it takes all positive, it 
is shown that it is possible to find a circular orbit even within the boson star. Details of the model are as 
follows: The mass for both, the black hole and the boson star, was taken as M = 0.633mp[/m, with m 
the mass of the constituent bosons, a free parameter. The central density, in usual dimensionless units, is 
equal to 0.27. No self-interaction was considered. All figures in this Section are based on this same model, 
but the behavior is generic. The radial coordinate is x = mr. 



relativistic non-baryonic potential, generated by a non-rotating compact object, there are circular orbits for 
every possible value of the radial coordinate, including those which are inside the structure. It is interesting 
to note, in addition, that < 1, so there are no unbound orbits, hyperbolic in energetics ||3^. 



The circular orbits might not all be stable. Stability requires that, when evaluated using the general 
expressions given in Eq (13), V{r)" < 0. This yields 



V{r)"{2B - rB') = -GBB'r^ + 4{B'yr^ - 2BB"r^ < 0. (14) 

Again, this reduces to the Schwarzschild case when B{r) = 1 — 2M/r. The results of this computation in a 
generic boson star potential are shown in the left panel of Figure 6. For all values of r, V{r)" < 0. Then, 
in a relativistic non-baryonic potential, generated by a non-rotating compact object, all circular orbits, even 
those within the structure, are stable]^ 



^An alternative derivation of the previous results can be obtained as follows. Consider again the equation of motion for r, 
written in the following slightly modified form: 



By multiplying both sides of the latter equation by grrgtt, and defining a new radial coordinate by z = J y/grrgttdr, we get 
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Figure 6: Left: Radial dependence of the second derivative of V{r), both for a Schwarzschild black hole 
(dash line) and a boson star (solid line). The minimum radius at which a stable orbit exists in the black 
hole case is r = 6M, it is marked with a crossing in the plot. For the boson star, all circular orbits are 
stable. Right: Effective potential for different values of L/M, with M the dimensionless total mass of 
the star, as in Figure 5. Two cases of the black hole effective potential are shown for comparison. BS(H) 
stands for boson star (black hole). The x-axis is the z variable, introduced in the text. For the case of a 
black hole, z and x coincide. 



5.2 From relativistic orbits to relativistic spectrum 

The power per unit area generated in a disc rotating around a compact central object, which produces a 



relativistic potential, is given by |37] 



Ml 

D{r) = 

47rr 



2M\ 

B{1 

r / 



dr 



1.4706 X l^-^'^m^GeV]^^^^^^ x 



cm^s 



B{x) 1 



47r 
dx 



- nL^) ' r [e^ - ni^) (^] dx. (17) 



= — 2Veff{r). We have considered that r = t{z)^ being the efTective potential defined by 



(16) 



This potential is depicted, for different values ol L/AI (with M being the total mass of the object in dimensionless units) in 
the right panel of Figure 6. Classical mechanics allows to extract, in the usual way, the orbital behavior. Superposed in the 
same plot of Figure 6, we have also depicted the black hole effective potential. Differences between them arise just from the 
rr-metric coefficient B(r), and are noticed in the innermost regions. A particle with any given energy can, in the non-rotating 
boson star vicinity, be in an stable, circular orbit, at any value of the radial coordinate. Depending on its energy, it can 
encounter one or two turning points, but there is no capture trajectory, consistently with the absence of singularity. 
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Dimensionless radial coordinate Dimensionless radial coordinate 

Figure 7: Left: Integrands for D(r), both for a black hole and a boson star of the same mass. The BH last 
stable orbit is marked. Right: Relativistic power per unit area. In dash lines, we show, for comparison, 
also the Newtonian results for both cases. 



All physical quantities are to be used in dimensionless form (e.g. B — > B{x)), and are explicitly given by 

E^j^^f, y\ (18) 

\2B-xB') ' \2{2B-xB')J ^ ' 

Clearly, when comparing with the output produced by a disc rotating upon a central black hole, D[r) 
will change not only because of the modifications in $7, E^^ , and L^, but also because of the change in 
the integration limits. Here, rms is the position of the innermost stable circular geodesic orbit. Using the 
latter expression and all obtained results for other physical parameters of the particle orbits, we can get 
the relativistic results we were searching. 

Figure shows the results of some crucial intermediate computation needed to get the relativistic 
spectrum, and, on the right panel, the power per unit area in the relativistic disc. It can be seen there that 
there is an interesting effect produced by boson stars as central objects: While a relativistic treatment of 
the disc reduces the expected emission from the innermost parts of a disc rotating around a Schwarzschild 
black hole, it dramatically enhances the expected one when a boson star is in the scenario. This has the 
consequence of a hardening in the emission spectrum, which is shown in Figure ^. The deviations shown 
between the spectra are impressive. Table 3 shows some of the values used to construct this latter figure. 



6 Conclusions 

In this paper we have modeled a very simple accretion disc rotating around a static supermassive boson 
star, although we have given the scaling property that shows how to extend these results to other mass 
domains (equivalently, to other single boson mass cases). The disc was assumed steady, with a constant 
accretion rate, and thin, so that the standard theory can be applied. Throughout the paper, we have 
made a comparison of all results with those obtained for discs rotating around Schwarzschild black holes of 
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Figure 8: Relativistic emission spectrum for an accretion disc rotating around a static boson star. The 
black hole case is shown for comparison. 



Table 1: Comparison of relativistic spectra. 



V [1/s] 


vL(v) [erg s ^\ 


uL{u) [erg s ^] 




BS 


BH 


IQl'^ 


2.00 X 10^5 


1.41 X 10^5 


10^2 


2.00 X 10^1 


1.41 X 10^1 


1014 


1.77 X 10^7 


1.18 X lO^'^ 


10^6 


1.36 X 10^^ 


3.46 X 10^3 



the same mass. Our aim was to see whether the emissivity properties of the accretion disc are noticeably 
changed when the central object is. More complicated models for the accretion process as well as more 
realistic models for the star (as those in which the star is rotating) can be considered. We hope this work 
will encourage further analysis. 
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